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StJMMARY 


This paper presents a simple solution to the problem of the acoustical 
coupling between a rectangular structure, its air content, and an external 
noise source. This solution is a mathematical expression for the normalized 
acoustic pressure Inside the structure. The paper also gives numerical results 
for the sound-pressure response for a specified set of parameters. 


INTRODUCTION 


The formulation of the problem is bas,ed on the following assumptions; 

1. The structure consists of a three-dimensional chamber, oriented with 
respect to a Cartesian coordinate system as shown in figure 1. The boundaries 
of the chamber are rigid except for an elastic wall, of homogeneous material, 
exposed to an external noise source and clamped at all edges . 


2. The external noise source is assumed to be a pure- tone |i.e., single- 
frequency] signal of known amplitude and frequency. 

3. The air inside the chamber is considered to behave as a compressible 
viscous fluid tindergoing oscillations of small magnitude. 

4. The elhstie wall is considered to behave as a vibrating plate with 
linear damping. 

The external^ incident noise-pressure disturbance causes the elastic wall 
to vibrate in the transversal direction. Inducing pressure fluctuations inside 
the chamber with a subsequent internal pressure loading on the elastic wall. 
The solution for the acoustic pressure inside the chamber, when damping and 
viscous effects are neglected, has been presented in reference 1; the effects 
of air viscosity and wall damping are included in the analysis given in this 
paper. 
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SYMBOLS 


a,b ,c 

^ij 

C 

D 

E 

f 

h 

1 

^mn 

P 

Pc 

Pc/2 

Po 

Po 

R 

t 

W 

w 

x,y,z 

a 


Y 

5 


height, width, and length of the chamber 

coefficients in the expression for the elastic-wall deflection 
speed of sound 

elastic-wall bending stiffness j D = Eh^/12(l-o^) 

Young's modulus for the elastic wall 

a function of a, 6 , y , and 6 , defined by equation (28) 
thickness of the elastic wall 

±2 =: _i 

air viscosity damping coefficient 
constants of integration 

coefficients in the expression for the acoustic pressure 
sound-pressure level inside the chanber 
sound-pressure level at z *= c 
sound-pressure level at z = c/2 

external sound-pressure level acting on the elastic wall 
amplitude of the time-harmonic^ external sound-pressure level 
weighting fionction used in the weigh ted- residual method 
time 

components of air velocity inside the chamber 

deflection of the elastic wall in the positive z-directlon 

mode shape of the elastic wall 

Cartesian coordinates 

wall-to-air mass ratio 

wall-to-alr stiffness ratio 

wall-to-air interaction damping ratio 

dimensionless air viscosity damping 


4 4 4 

^ + 2 ^ ^ 
3x^ Bx^ay^ 3/ 


T1 

X jP >V 


wall damping coefficient 
chamber width-to-h eight ratio 
separation constants 

Separation constants in the expression for the acoustic pressure 
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5 

p 

0 


dimensionless parameter defined by equation (29-b) 
density 

Poisson^ s ratio for the elastic wall 
0 ) circular frequency of the external noise 

Subscripts: 

a refers to the air inside the chamber 

max denotes "maximum value" 

w refers to the elastic wall 

Superscript: 

- .refers to dimensionless quantities 


MATHEMTICAL FORMULATION 


The governing dynamic equation for the elastic wall is: 

n4tr ^ ^ 3W h 9^W _ 1 , . 

^ ^ D 9t‘^ ^WD 2 D 

dt 

The acoustic wave equation for the air contained in the chamber is: 

^ + 1 ^ 2 . = —L. ( 1 ^ + Ij. 

2 ^ 2 ^ 2 2 ^ 2 ^ a 9t“^ 

9x 9y 9z Ca 3t 

The boundary conditions for the problem are as follows: 

a) The edges of the elastic wall are considered to be clamped: 
W(0,y,t) = W(a,y,t) = W(x,0,t) = W(x,b,t) = 0 


( 1 ) 


( 2 ) 


(3) 


W 

9x 


(o,y,t) 


9x 


(a,y,t) 


f A r f\ n 

^(x,0,t) = -^(x.b.t) ^ 0 


(4) 


b) The normal component of the internal air velocity near a rigid boundary 
is zero: 

Vjj ( 0 ,y , z , t) = (a ,y , z , t ) = (x , 0 , z , t) =v (x ,b , z , t) =v^ (x ,y , 0 , t) =0 (5 ) 

c) The normal component of the internal air velocity near the elastic wall 
is equal to the wall velocity: 
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V2(x,y,c,t)^ (6) 

The relationship between the internal air pressure and components of internal 
air velocity are: 


9v 1 r, 9v T 

--i=- JL^.Wv — 


Pg^ 9x ax 8^ 


Pa ay ’ 9t 


- k v„ (7) 


p 3z a z 
a 


The external noise pressure is considered to be harmonic in time and expressed 
by : 

•o /ON 

p = P e (8) 

*^o o . 

and the objective is to find: 

W = W(x,y,t) and p p(x,y,z,t) (9) 

The solution of equation (2), by separation-of- variables technique, is: 

00 oo 

p(x,y,z,t) = ^ k Cos-2^ Cos^-^ Cos v z (10) 

^0 '^0 a b mn 

2 

„ w -imk 2 2 

where - <— ) - ( 5 -) ( 11 ) 

o 

a 

Application of equations (6) and (7) yields; 


■|H . . — r r K V CoeSS CosSS- Sta V c 

3t p (k +xoi) ^ ^ mn mn a b mn 

a a niFO n=0 


SiT 7 

The value of -rr is found by integrating equation (1) under the loading 

O t 

conditions as indicated below: 


i M + k i-W 

D 3t '■^wD 3^2 


itotn w ® 

==; — Cos®^ Cos^^ Cos v„c) - P- 


For a solution of the form: 


W(x,y,t) = w(x,y)e 


equation (13) reduces to: 

,i 2 h 1 

V® + — W- jw = J 


S» = Frf5irK CosSS CosSa Cos v^c\ - P„1 (15) 

D D \^— mn a b mn y o 

*- 111=0 n= 0 ' -* 
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Galerkin’s method is used to find an approximate solution to equation (15). 
In this method, an approximate solution which satisfies the boundary con- 
ditions of equations (3) and (4) is first assumed as follows; 


M' N 


1 . 


[1 - 1 

a 


b 


(16) 


Coefficients a±j are found such that equation (16) satisfies equation (15) 
and the sum of the weighted residuals is identically zero over the region 
of integration, i.e.; 


b 


^ 0 


, io)C - , , r 00 00 1 

7 W + - Jw- p w - -(j^, £^KffinCosS^ Cos^^ Cos 


R dxd3F=0 

(17) , 


This process is known as the weighted residual method and R is the weighting 
function. In Galerkin's method, the weighting function is made eqtial to the 
shape fxmction defining the approximation. In general, this leads to the best 
approximation when 


R = 

- Coa'2i-«2**l 




a 


D 


(18) 


Equation (17) applies for every pair of- integers i aiid j . Generally there 
are M*N simultaneous equations of this form to be solved for the coefficients 
a^j . Por the case of low-frejquehcy normally- incident external noise only the 
diaphragm motion of the wall (first mode) will be excited. For this case, 

M = 1 and N = 1 , and integration of equation (17) leads to: 


^22 I ^20 V 

— Cos V 22 C - — Cos V2qC - -j- Cos Vq2C + Cos v^qC 


- P 


and W(x,y,t) = a^^ (1 - Co s ~ -^ ' ) (1 - Co s^^) e^^^ 


(19) 


( 20 ) 


The values of ^22*^20*^02* %0 found by substituting the deflection 

from equation (20) into equation (12): 


(-“ 


+ ik w) p a , 1 
a a 11 


(l-Co,?S) (l-c„s^) =E I K, 


m=0 n=0 


V Cos'^^^ Cos-^ 
mn mn a b 


Sin V c 
mn 


( 21 ) 


When the left and right sides of equation (21) are equated term by term, it 
is found that all the constants are zero except Kqo,Kq2,K2o> ^22* 
These are easily found and substituted in equation (19) to give: 
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(22) 


4v 


22 


2v 


20 


2v 


02 
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00 


Referring to equation (10), the acoustic pressure inside the chamber can now 
be written as; 


p(x,y,z,t) = e'"*”^ ^Kqq Cos VqqZ + Kq 2 Cos^^ Cos 

K„„ Cos^^ Cos + K„„ Cos'^^ Cos'^^ Cos v„„z 

20 a 20 22 a b 22 


(23) 


To generalize the solution obtained above, the following dimensionless 
quantities are introduced: 


- X - y 
X = — , y = > 
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r _ _w n - -S 
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(24) 

(25) 

(26) 


Using the above dimensionless quantities , equation (23) can be written as: 


icut 


p(x,y,z,t) = e f(a,3,Y,5) 
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(28) 


and 

a = p h = 


k c 



a 


fl _ (1 + 1.056n^)pA^ - ^ S_ 

» p _4 > Y f, fi 

0.0284b a 

= k = i r = 0.9221(3n^ + 2n^ + 3) 

^ ^ ’ (1-0^) (1 + 1.056n^) 


(29-a) 

(29-b) 


Equations (2^, (28), and (29) constitute the analytical solution to the acousto- 
structural problem. These eqtiations show that the normalized pressure dis- 
tribution within the chainber is a harmonic function of time and depends on the 
following dimensionless parameters; 

1) wall- to- air mass ratio, a 

2) wall-to-air stiffness ratio, 8 

3) wall-to-air interaction damping ratio, y 

4) dimensionless air viscous damping, ^ 

5) dimensionless frequency, u) 

6) dimensionless space coordinates, x,y,z 

7) enclosure dimensions, a,b,c 


NUMERICAL RESULTS 


To obtain quantitative values of sound-pressure level as a function of exteraal 
noise frequency, a cubical chamber (a=b=c) is assumed, and the amplitude of 
the dimensionless sound pressure at the center of the chamber (x=y=z=c/2) is 
found: 


^c/2 

max 


Numerator 

Denominator 


(30) 
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The response of this "advanced" three-dimensional model, as given by equation 
(30) , is compared with that of a "simplified" one- dimensional model obtained 
by replacing the elastic wall by a simple spring-mass system. For this sim^ 
plified model the amplitude of the sound pressure at the center of the cubical 
chamber is ; 


1 



If the effects of wall damping and air viscosity are neglected, the results 
given by equations (30) and (31) agree with the solution in reference 1, in 
which damping effects were not considered. Figures 2 and 3 show the frequency 
response, for a particular set of dimensionless parameters a and 3 , over 
the audio-frequency range, for the special case of y “ 0 snd 6=0, i.e. , 
when damping effects are neglected. These figures show that at intermediate 
frequencies and at the high-frequency end of the audible spectrum, the pre- 
dictions of "advanced" and "simplified" models are quite similar. 

When damping effects are included, i.e., when both y and 6 are not zero, 
the digital computer program for the frequency response is very complicated, 
involving complex numbers and requiring double-precision (16 digits) accuracy. 
Results for this case will be published later. 
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Figure 1.- Three dimensional model of sound transmission through an 
elastic wall to a rectangular chamber. 
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DIMENSIONLESS PRESSURE,(p, K) DIMENSIONLESS PRESSURE, (p„). 





